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Abstract 

In this note, we present the formalism to start a quantum analysis 
for the recent billiard representation introduced by Damour, Henneaux 
and Nicolai in the study of the cosmological singularity. In particular 
we use the theory of Maass automorphic forms and recent mathemat- 
ical results about arithmetical dynamical systems. The predictions of 
the billiard model give precise automorphic properties for the wave 
function (Maass-Hecke eigenform), the asymptotic number of quan- 
tum states (Selberg asymptotics for PSL(2,Z)), the distribution for 
the level spacing statistics (the Poissonian one) and the absence of 
scarred states. The most interesting implication of this model is per- 
haps that the discrete spectrum is fully embedded in the continuous 
one. 
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1 Introduction 



The beautiful singularity theorems of S. Hawking and R. Penrose [TT] es- 
tablish the presence of a singularity in the generic solution to Einstein field 
equations under very reasonable assumptions of the kind of matter. These 
theorems use tools from differential topology (Morse theory for Lorentzian 
manifolds) but do not say anything about the behavior of the singularity, 
being existence results. Thus the works by V.A. Belinskii, I.M. Khalatnikov, 
E.M. Lifshitz (BKL in the following) [2] regarding the general behavior of a 
homogeneous cosmological singularity are very important in understanding 
the (space-like) singularity in general relativity. Their approach is based on 
the assumption that close to the cosmological singularity the spatial gradi- 
ents can be neglected with respect to the temporal ones in the field equations, 
this being equivalent to a spatial decoupling of the points on the singular 
hyper-surface (BKL limit). Thus each point (i.e. the space part of the met- 
ric) evolves according to an infinite succession of Kasner eras, leading to 
an approximate discretized dynamics which turns out to be isomorphic to 
the Gauss map (the continued fraction expansion of a real number), hence 
chaotic. The analysis of BKL has been recently improved and generalized 
by T. Damour, M. Henneaux and H. Nicolai (DHN in the following) [6] 
to the case of a general inhomogeneous cosmological singularity and super- 
gravity/string theories (extra dimensions, exotic matter etc). Their result is 
that, in the BKL limit, the dynamics of pure gravity in 4 dimensions can be 
reformulated as a billiard motion in a suitable domain embedded in a flat 
3-dimensional pseudo-Riemannian manifold (we comment on the chaoticity 
of this billiard too). This billiard table has a basis on the hyperbolic plane 
corresponding to the fundamental domain of the extended modular group 
PGL(2,Z), which is a so-called arithmetic group. The arithmetic nature of 
the asymptotic billiard allows to derive a precise quantum analysis and make 
some predictions for this model using powerful and recent theorems about 
the quantum behavior of arithmetical dynamical systems. This quantum 
analysis is the main point of this note. In particular, we can give the precise 
form of the wave function, count the asymptotic number of quantum states 
and comment on the so-called scarred states. Finally, the role of the Selberg 
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trace formula for PSL(2,Z) as a convergent semi-classical quantization rule 
for gravity in this regime is briefly pointed out. 

Two appendices, on the Kac-Moody algebra HA^ and the theory of Maass 
waveforms, are added at the end and should be read in parallel with the 
relevant sections. 

This paper is based on [TJ, which also contains more details, references and 
background material. For the BKL limit, the original references and many 
other related questions we suggest the recent reviews [18j and [13]; the latter 
deals especially with the cosmological billiards in relation to other theories 
(supergravity etc) too. 

We want to stress that, usually, with the term chaotic cosmology one refers 
to the fact the asymptotic (i.e. in the limit towards the singularity) classi- 
cal evolution exhibits some chaotic features, but one has to consider what 
are the implications of that in the quantum behavior of the system. As we 
will see the corresponding quantum system is not "chaotic" , since the quan- 
tum manifestations of classical chaos for arithmetical systems are closer to a 
classically integrable system, and this is due to the arithmetic nature of the 
billiard. In fact, the existence of Hecke operators, which form at the quantum 
level an infinite family of commuting self-adjoint operators which mimic an 
integrable system, allows for an anomalous spectral statistics (a Poissonian 
law for the high energy spectrum of the quantum Hamiltonian instead of a 
GOE/GUE ensemble) and allows to derive the quantum unique ergodicity 
theorem. 

2 The BKL oscillatory behavior and the bil- 
liard representation 

Let us briefly review the BKL result, restricting to the case of a homogeneous 
Bianchi IX universe, or Mixmaster as renamed by C. Misner for the chaotic 
features it exhibits. Moreover, we consider only the case of pure gravity in 
4 dimensions. Then, the spacetime metric can be written in a synchronous 
reference system as 

ds 2 = -dt 2 + dl 2 = -dt 2 + h ij dx i dx j (1) 
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with hij diagonal 

hij = a 2 (t)lilj + b 2 (t)mimj + c 2 (t)ni7ij (2) 

and where the frame vectors l,m,n depend on the spatial coordinates x l . 
The behavior of the three scale factors a, b, c has been studied by BKL in the 
limit t — > C0. Their evolution is represented in the picture. The asymptotic 




Figure 1: Evolution of the logarithms of the scale factors (a = In a, /3 = 
In b, 7 = In c) in terms of a logarithmic time Q = — In t for a Bianchi IX 
universe (Mixmaster). 

dynamics is such that there exist an infinite number of Kasner eras, each one 
formed by a certain number of Kasner epochs. Precisely, during each Kasner 
era, two of the scale factors oscillate and the third one decreases monoton- 
ically. On passing from one era to the next one, the monotonic decrease is 
transferred to another of the three scale factors. The time interval in which 
each scale factor is represented as a straight segment is called a Kasner epoch, 
i.e. it is a regime of the solution to the field equation where one can approx- 
imate the exact solution with a Kasner solution. Given an irrational number 

1 Regarding the validity of the BKL regime, we have to say that it should be settled 
down in the pre-inflationary Universe, thus unaccessible to present observation. Thus 
it concerns the very early stages of the Universe evolution, matching the Planckian era 
with the inflationary behavior. Of course, one should always consider the problem of the 
falsification of a theory of quantum gravity, since it deals with Planck scale physics. 
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u > 1, the number of Kasner epochs contained in a Kasner era corresponds 
to the partial denominators in the continued fraction expansion for u. 
In other words, the result is that the asymptotic evolution can be described 
as infinite succession of Kasner epochs with a certain law of replacement of 
the Kasner exponents when passing over from one epoch to the next one. 
This law corresponds to the Gauss map 

(3) 
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whose metric entropy can be computed 
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h(T Gauss ) = ^-^ (4) 

Note that the Gauss map accounts only for the transitions between successive 
Kasner eras (i.e. in the picture between the times Qq and Q\ and so on). One 
can do better and consider also the oscillations in two of the scale factors 
inside a single Kasner era and include them in a more complete map which 
describes the discrete evolution (still approximated). This can be realized 
through the Farey map, whose topological entropy is given by 

T(T Farey ) = 2 In 2 (5) 

This map accounts for oscillations (one pair of axes oscillates while the third 
one decreases monotonically) and bounces (when the role of the three axes 
are interchanged and a different axis decreases monotonically) , according to a 
chaotic Farey tale j5]. This paper uses fractal techniques (which are observer 
independent) to claim that the Mixmaster universe is indeed chaotic. 
Finally, note that in the BKL approach (although fascinating) it is not easy 
to write down explicitly a quantization procedure. 

BKL also considered the case of a inhomogeneous universe and showed that 
there exist a generalized Kasner solution near the singularity 

dl 2 = dx % dx^ (6) 

= a 2 lilj + b 2 mimj + c 2 niUj (7) 

with 

a~t Pl , b~t Pm } c~t Pn (8) 



4 



and 

Pl(j) + Pmtx 1 ) + Pn (x l ) =p 2 l (x i ) + + pKx') = 1 (9) 

but now, differently from the homogeneous case, the Kasner exponents are 
functions of the spatial coordinates. One expresses this fact by saying that 
inhomogeneous cosmology can be modelled on a Bianchi IX universe, since 
there exists an analog infinite succession of Kasner epochs/eras with local 
behavior, i.e. each point on the singular hypersurface evolves with the same 
laws (Bianchi IX) but with different parameters. In this sense, one says that 
the behavior of a generic singularity is local and oscillatory. 
The positivity of the metric entropy for the Gauss map implies its ergodicity. 
This result was derived for the first time by Artin (much before the introduc- 
tion of entropy in dynamical systems) relating the properties of the Gauss 
map to the geodesic flow on PSL(2,Z)\H, which was known to be ergodic 
(see [21J). Thus, one can ask if there exists a precise physical system associ- 
ated to this geodesic flow in the same way as the asymptotic evolution of a 
Bianchi IX universe can be approximately described by the Gauss map. The 
answer to this question is positive and the physical system is pure gravity in 
4 dimensions near the singularity and without any symmetry assumption for 
the metric. This leads us to the billiard representation introduced by DHN. 
The result by DHN [6] is that the BKL oscillatory behavior of pure grav- 
ity in 4 dimensions can be reformulated as a billiard motion in an auxiliary 
spacetime^l More precisely, the dynamics of the generic solution to Ein- 
stein equations close to the cosmo logical singularity (and in the BKL limit, 
i.e. assuming BKL's conjecture) is equivalent to a null geodesic motion 
(i.e. on a light ray) inside a billiard table given by a Coxeter polytope in 

2 One can speculate whether this representation is photographic or holographic (since 
we go from a 4-dimensional spacetime to a 3-dimensional auxiliary Minkowski space), 
according to the words of T. Damour at the last 11th Marcel Grossmann meeting. 
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a 3- dimensional Minkowski spacej The remarkable fact is that this poly- 
tope corresponds to the Weyl chamber of the hyperbolic Kac-Moody algebra 
HA^, the canonical hyperbolic extension of A x (su(2)). The walls are the 
hyperplanes pointwise fixed by Weyl reflections in the three simple roots 
which define the algebra and the reflections at the walls are elastic. In this 
language, each Kasner epoch is represented by a null geodesic segment be- 
tween two successive reflections. In particular, given a Kasner epoch and the 
wall where this epoch crashes/ends, the following one is obtained by Weyl 
reflection with respect to the simple root orthogonal to that face of the Weyl 
chamber. In other words, Weyl reflections with respect to simple roots send 
null geodesic segments into null geodesic segments, i.e. transform Kasner 
solutions into Kasner solutions (with different values of Kasner exponents). 
Because of this, we have a hint that the Weyl group of the algebra acts on 
the space of solutions of Einstein equations, transforming Kasner solutions 
into Kasner solutions which can then be associated to the simple roots (and 
more generally to the real roots of the algebra). Of course, this would be 

3 Generic means without any symmetry assumption for the metric. If accidental symme- 
tries are present in the metric, the analysis is not valid. For example, it is known that the 
Schwarzschild solution 

ds 2 = - (l-^dt 2 + (l-^j 1 dr 2 +r 2 dn 2 (10) 

has a spacelike singularity at r = 0; moreover, inside the horizon (r < 2m) the r coordinate 
is time- like (there is a minus sign in front of dr 2 ). If we take the limit r — > 0, we obtain 

lim ds 2 = — dt 2 - — dr 2 + r 2 dQ 2 (11) 

r^o r 2m 

which is a Kasner metric 

- dr 2 + T - 2/3 da 2 + T 4/3 (sin(6>)#) 2 (12) 

once we put r = |^j=, a = (4m/3) 1/3 t, 6 = (9m/2) 1 / 3 6> and = (9m/2) 1 / 3 0, i.e. the 
Schwarzschild solution corresponds (in the neighbourhood of the singularity) to a single 
Kasner epoch, not to a never-ending succession of Kasner eras. 

Finally, the analysis does not apply to time-like (like the one in the Reissner-Nordstrom 
solution, the charged black hole) or null singularities where a causal decoupling of spa- 
tial points does not occur. The question about the general behavior of not-space-like 
singularities is still open. 
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a technique to generate solutions (also approximate) to Einstein equations 
starting from very simple solutions. This was already done by R. Geroch, 
who found that gravity reduced to 1+1 dimensions has a hidden A^ 1 sym- 
metry (affine su(2)), which is a sub-algebra of ELA| . Thus in this case, 
the affine algebra (or better its corresponding infinite-dimensional Lie group, 
known as the Geroch group) can be considered a symmetry of the theory. 
One can go further and check if the whole hyperbolic algebra (not only its 
Weyl group, on which we focus in this work) is a symmetry of gravity and 
in which regime, see |13j . 

In this formalism only null geodesies are physical and correspond to Kasner 
solutions. Space-like and time-like geodesies seem to play no role. Note that 
the walls of the Weyl chamber are time-like, that is their orthogonal vectors 
are space-like, since for the simple roots one has (aij, c^) = 2 > 0. Because 
of this, every reflection preserves the null character of the velocity vector. 
The (flat) metric in this Minkowski space is given by the Cartan matrix of 
the algebra 



whose signature is ( — h+). For this reason, one has three kinds of geodesies, 
but only null geodesies appear. 

Note that in the DHN approach, the 3-dimensional billiard is peculiar, in 
fact the walls move together with the trajectories. During this motion the 
geometry of the billiard remains the same, i.e. the structure of the algebra 
is preserved and the Weyl group is uniquely determined. This means that 
the walls of the (reduced) billiard on the hyperbolic plane are motionless and 
the 2-dimensional analysis is much easier than studying the 3-dimensional 
billiard with (say) co-moving walls. As we will say many times in this work, 
the Weyl group of the algebra is an arithmetic group (being isomorphic to 
PGL(2,Z)), thus we are coping with a physical system represented by an 
arithmetical dynamical system. 

The very essence of arithmetical chaos (see [I]) must be researched in its 
quantum aspects, since from the classical point of view no difference exists 
between generic chaotic dynamical systems and arithmetical dynamical sys- 




(13) 
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terns, except for the degeneracy of lengths of periodic orbits. In fact, for 
generic systems, one does not expect a degeneracy of this kind but the ones 
which come from the symmetry of the model. For example, for time- invariant 
systems, this multiplicity, on average, is 2, which corresponds to the two pos- 
sible ways to run a geodesic. Arithmetical systems are very exceptional, since 
the degeneracy of lengths of periodic orbits grows exponentially 

1 e L ' 2 

cm (14) 

where L is the length of such a geodesic and C is a constant depending 
of the particular group (C = 1 for PSL(2, Z)). To this, one should add 
the Horowitz-Randol theorem which states that these degeneracies are un- 
bounded. The large degeneracy of lengths of periodic orbits seems to have no 
importance in classical mechanics. These systems are ergodic as any other 
model on negatively curved surfaces, but the quantum aspects are anomalous. 
Thus, it is very interesting to start a quantum analysis of the cosmological 
billiard using its arithmetic nature. Before doing that, in the next section 
we comment on the (classical) chaoticity of the cosmological billiard. 



3 General Considerations on the Cosmologi- 
cal Billiard 

Let us stress again that the billiard motion occurs in a pseudo-Riemannian 3- 
dimensional spaced, not in a Euclidean space. This space is f}Jj, in the notion 
of Kac [Mj (see also Appendix [A]) . Thus the billiard motion is an interrupted 
null flow on a pseudo-Riemannian manifold. This situation is different from 
the the typical billiards which are embedded in Riemannian manifolds. In 
fact, in pseudo-Riemannian manifolds, one has three kinds of geodesies, and 
consequently three kinds of geodesic flows (indeed one first needs to define in 
a rigorous way flows for manifolds with an indefinite metric). The second step 

4 We prefer the term pseudo-Riemannian instead of Lorentzian because for people who 
study billiards as dynamical systems, the word Lorentzian has a completely different mean- 
ing. 
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is the definition of the billiard AoWlI, which is even more involved. In fact, 
because of the three kinds of directions (null-like, time-like, space-like), each 
reflection on a wall depends both on the character of the incoming trajectory 
and on the character of the wall. 



Figure 2: Schematic representation for the cosmological billiard in the Weyl 
chamber of HA^ and projection on the Poincare disc. 

To clarify this important point, let us consider the usual billiard motion in 
a convex polyhedron Q in a Euclidean space, that is a Riemannian manifold 
with the usual (diagonal) flat metric (see [2T] ). The boundary of this poly- 
hedron is given by the union of its faces IV Each elastic reflection in the wall 
Tj is given by the geometric reflection <jj 



where v is the velocity vector of the incoming trajectory, rij the unit vector 
orthogonal to each face and (, ) is the standard Euclidean scalar product 
(compare with the expression of the Weyl reflections in Appendix |Aj). To 
understand the properties of the billiard flow inside Q, one has to consider 
the group Gq generated by the reflections this is a subgroup of all 

isometries of the Riemannian manifold. The ergodicity of the billiard depends 

5 The problem is the definition of the billiard law, i.e. how to use the rule "angle of 
incidence = angle of reflection" in a pseudo-Riemannian setting. When one speaks of 
billiards, one must be careful, because together with the billiard table, one has also to 
specify the billiard law. In Euclidean spaces, one usually considers the law of geometrical 
optics, but in principle one can study different billiard laws (in the same table), as it 
happens for outer billiards which have nowadays surprising applications. 





(Ti{v) = v - 2(ni,v)ni 



(15) 
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on Gq, precisely: if Gq is a finite group, then the billiard flow inside Q is not 
ergodic (thus it is not chaotic). In 2 dimensions, the finiteness of the group 
is equivalent to the commensur ability of all angles on the polygon Q. 
The situation for generic polyhedra is still open. Indeed, one knows that 
the metric entropy of a billiard inside an arbitrary, not necessarily convex, 
polyhedron is zero. All these statements are valid for Riemannian manifolds. 
In our case, the analogy is the following. The polyhedron Q is replaced by the 
Weyl chamber of HA^, the Riemannian manifold is replaced by the pseudo- 
Riemannian space fj^ and the billiard motion occurs on null segments. The 
faces of the billiard are timelike walls. Finally, the reflection group Gq is 
replaced by the Weyl group of the algebra W, which is an infinite Coxeter 
group. At the moment of writing this work (and modulo over-present igno- 
rance), there are no mathematical results which imply the ergodicity (or the 
integrability) of our billiard motion, because the difficulty is that the motion 
occurs in pseudo-Riemannian manifolds, thus we can not invoke the previous 
theorems for billiard flows in Riemannian manifolds. In particular, we know 
that W is isomorphic to PGL(2,Z), and that the free motion on the hy- 
perbolic plane inside a region with finite (hyperbolic) area is a chaotic flow, 
being an Anosov flow (|21j). Thus we can state that there exists a chaotic 
motion on the "basis" of the billiard, but this does not imply the chaoticity 
of the null billiard flow in the full 3-dimensional Weyl chamber. In fact, the 
full motion could be less chaotic or even integrable. For example, the Buni- 
movich stadium is a well-known example of a 2-dimensional Euclidean chaotic 
billiard, but if we consider a 3-dimensional Euclidean billiard raising the sta- 
dium as a basis in the z— direction, then the billiard is integrable, because of 
a translation symmetry in the z— direction. Thus, a billiard with a chaotic 
basis (i.e. which has the property of being ergodic or mixing or hyperbolic 
etc) is not necessarily chaotic. Our case is even more complicated, because 
the full 3-dimensional billiard is a Coxeter polytope in a pseudo-Riemannian 
space, and the projected billiard corresponds to the fundamental domain of 
PGL(2, Z) on the hyperbolic plane (which is a Riemannian manifold) §. 



6 With the projection of the 3-dimensional billiard on the hyperbolic plane, we mean 
that we consider the geodesic segment on the hyperbolic plane corresponding to the inter- 
section in 3 dimensions of a null segment with a plane through the origin. 
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To conclude this section, one can not state that the Kac-Moody billiard is 
chaotic only from the property that the billiard table is contained inside the 
light-cone, but more powerful (and at the moment missing) theorems are 
needed. The study of billiard flows on pseudo-Riemannian manifolds is a 
completely new field and should be an interesting topic of future research. 
It seems that the only mathematical papers on this subject are very recent 
and due to B. Khesin and S. Tabachnikov [16], where they define a billiard 
law which is area-preserving and does not change the type of a geodesic. 
Let us also observe that the metric entropy of the geodesic flow inside the 
fundamental domain of PGL(2, Z) can be explicitly calculated^ 



thus the billiard representation on the hyperbolic plane and the BKL ap- 
proach based on the Gauss map are not equivalent, as they are not isomorphic 
as dynamical systems. In fact, the first describes the behavior of a generic 
inhomogeneous singularity, the second the fate of a Bianchi IX homogeneous 
universe. Yet, as we mentioned above, a famous result due to Artin relates 
the ergodicity of the geodesic flow on PGL(2, Z)\H to the ergodicity of the 
Gauss map. Artin's theorem supports the conjecture that the behavior of 
a generic space-like singularity is somehow well described by a Bianchi IX 
homogeneous cosmological model, conjecture for which nowadays there is a 
lot of numerical evidence [3J. At this point, we should add that some the- 
oretical progress has been made in understanding the general behavior of 
the cosmological singularity thanks to C. Uggla et al (see [22] for a review), 
who have shown that asymptotic silence holds, i.e. particle horizons along 
all timelines shrink to zero for generic solutions. With Uggla's words at the 
last 11th Marcel Grossmann Meeting, "Everybody dies alone". 
We can make an interesting consideration when we pass from the 3-dimensional 
Minkowskian billiard to the 2-dimensional Euclidean one. The light rays 
which move in the auxiliary Minkowskian space have no mass of course. 

7 In fact, for the geodesic flow {S 1 *} on a surface of negative constant curvature K, the 
following formula holds 



h({S 1 } on PGL(2,Z)\M) = 1 



(17) 




(16) 



and one usually puts the Gaussian curvature JC 



1. 
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When we study the wave equation on the reduced domain (i.e. the Schrodinger 
equation for PGL(2, Z)), there is a subtlety. Indeed, this equation allows for 
a mass different from zero because it corresponds to the quantization of a 
classically chaotic particle on the hyperbolic plane (we always put m = 1, 
h = 1). But this particle "arises" from light-like rays in 3 dimensions. Thus, 
we have an example of a mechanism to generate mass on a Euclidean mani- 
fold from a null motion on a Minkowskian manifold. 

4 The Wave Function of the Universe 

Let us now come the quantization of the cosmological billiard. In the lan- 
guage of [6] , we can decompose the motion in the 3-dimensional billiard in an 
angular part inside a fundamental domain for PGL(2,Z) plus a radial part. 
The latter gives a temporal contribution to the full wave function, because 
the vertical direction in the Weyl chamber is a temporal one. Thus, with 
the term wave function in this section we mean the space part of the full 
wave functional which is solution of the Schrodinger equation on the hyper- 
bolic plane. How to arrive at this equation after some changes of variables is 
already described in many works (see [18J and references therein), we briefly 
remind the main steps below. But let us stress that the wave equation on 
the hyperbolic plane derives from the Wheeler-DeWitt equation for the 3- 
dimensional problem. We believe that the time-dependence of the solution to 
the full wave equation does not change the physical spectrum, which is thus 
contained in the 2-dimensional billiard. In [TH] (following previous works 
by Misner), the authors transform the domain for Bianchi IX using an ap- 
proximation which makes the domain integrable (they modify the lower arc 
of the billiard with a horizontal segment, but note that the latter is not a 
geodesic segment on the hyperbolic plane) and with the same area. The 
same approach was used for Bianchi IX in [8] , where the authors identify the 

8 To be more precise, one should also say that the full wave function of the system is, 
heuristically, the product of the wave functions in each space point on the hyper-surface. 
But since the classical behavior is local, we need only to study the quantum behavior in a 
generic point, where we approximate the physical evolution by the billiard representation. 
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fundamental domain on the hyperbolic plane specifiqj for this situation and 
try to resolve the corresponding wave equation. 




Figure 3: The domain on the hyperbolic plane corresponding to the asymp- 
totic dynamics for a Bianchi IX universe (from [8]). 

Here we keep the exact domain. The novelty with respect to previous works 
is our widespread use of the arithmetic properties of the billiard, which, 
although mentioned in some papers by DHN, have not been taken in consid- 
erations before, especially at the quantum level where they play a decisive 
role. 

To be more specific, using the notation of [8], the 3-dimensional quantum 
equation reads as follows: 

- jit^ ^( ln "» 2 >f + 7ihr^ B - u( "^ = (18) 

where A LB = ^ ^(sinhC ^) + — U is a potential term 

f/(p,C,0) = (lnp) 2 [p( 4cosh ^(^-l)] (19) 

9 They restrict to Bianchi IX without rotation of axes. If one considers also the rotation 
of axes, the fundamental domain is the one for To (2) [17], which is a subgroup of PSL(2, Z). 
It is reasonable for the generic/inhomogeneous case to expect a generic discrete group, and 
in fact PSL(2,Z) occurs. 
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p is the radial (time) variable, ( and <fi are the angular (space) variables. This 
equation is exact, in the sense that it is derived from the Bianchi IX metric 
without any other assumption^. 

If we know take the limit p — ► + (which corresponds to the case Q — > +00 
in the BKL variables), then the potential term vanishes inside and is plus 
infinite outside the triangular domain pictured in the figure. This makes it 
possible to factorize the solution as follows 

* = $Gc#(C,0) (21) 

and we consider only the equation 

- A LB iPi = (22) 

which after some redefinitions of variables corresponds to the spectral prob- 
lem for the hyperbolic Laplacian (see [H] and [8], in this latter the approx- 
imate temporal evolution for the 3-dimensional wave function is computed 
too). 

To go beyond, as we have seen, the general inhomogeneous case is modelled on 
a billiard problem inside the fundamental Weyl chamber of the hyperbolic 
algebra HA^ whose projection on the hyperbolic plane is PGL(2, Z) and 
not PSL(2, Z), thus we have to consider the spectral problem with Neumann 
boundary conditions 

-Aip = Eip 
V e L 2 {F 3 ,p) 

d n ipU 3 = (23) 

where JF 3 is the halved standard modular domain, A is the hyperbolic Lapla- 
cian and /i is the usual measure on the hyperbolic plane (the notation of 
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In particular, in the original Misner variables the Bianchi IX metric (TTJ) can be written 



MS 



ds 2 = -dt 2 + S 2 (i)(e 2/ %^V (20) 

where f3 =diag(/3 + + >/3/3— , /3+ — — 2/3+) is 3 by 3 traceless matrix and the a^'s are 

a basis of one-forms for SO (3). In terms of these variables, the term V in the potential U 
is V — -|Tr(l — 2e~ 2/3 + e 4/3 ). Finally, the variables /3±,o are transformed to hyperbolic 
coordinates p, £, (j> [5]. 
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this section is explained in Appendix [B] with more details). Thus solutions 
to the Neumann problem for PGL(2, Z) are given by even Maass cusp forms 
for PSL(2,Z). In fact, the spectrum of A on PSL(2,Z) is purely discrete 
on the space of odd functions, but on the even space there is a continuous 
spectrum given by the interval [7,00) and, what is less obvious, a discrete 
spectrum. Thus the space of Maass waveforms (which for PSL(2, Z) are also 
cusp forms, i.e. they are zero at the cusps) splits into odd/even forms under 
the symmetry R\ 

1>(-z) = ±1>(z) (24) 

Since the translation T(z) = z+1 belongs to the modular group, an even/odd 
Maass form can be written as ip(x) = ^2 n W n (y)e 2nmx . Considering the 
eigenfunction condition Aip + Q + t 2 ) ip = and the square-integrability 
condition, one arrives at a more explicit form for the Maass cusp forms for 
PSL(2,Z), that is 

<if)(z) = a^n)y l l 2 K u {2<n\n\y)e 2mnx (25) 

where K it {2n\n\y) are modified Bessel functions such that K it (y) << e~ 2wy 
for y — > 00. The coefficients a^,{n) are the Fourier coefficients of the ip{z). For 
even forms, a^(— n) = a^(n), while for odd ones a^(— n) = —a^(n). Finally, 
since for PSL(2,Z) (and for all so-called Hecke triangle groups) there is an 
obvious symmetry with respect to the imaginary axis, one can write 

*( I ) = |^jr,( a ™,)j»« (26, 

depending on whether ip is even or odd. The precise statement, then, is that 
the wave function of the universe is an even Maass cusp form for the modular 
group PSL(2,Z). Apart from the trivial eigenvalue E = and the corre- 
sponding constant eigenfunction ip , no other eigenvalues/eigenfucntions are 
known analytically. Some eigenvalues, calculated via numerical investiga- 
tions, are given at the end of Section [U 

Actually, more is true. In fact, since PSL(2, Z) is arithmetic and the cuspidal 
spectrum is very likely simple, one can diagonalize the hyperbolic Laplacian 
and the Hecke operators simultaneously. Thus it turns out that the wave 
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function is a Maass-Hecke eigenform. This is an even more interesting state- 
ment. In fact, Hecke eigenvalues are multiplicative, i.e. \(mn) = X(m)X(n), 
and this should put some conditions on the physical interpretation of them. 
In recent years, in string theory, a lot of work has been done on counting the 
entropy of black holes (see for example [T9.j for a review). It turns out that in 
some cases this entropy is counted by the Fourier coefficients of certain au- 
tomorphic functions. Thus, it is remarkable that the quantization of gravity 
in the BKL limit leads to a wave function which has automorphic properties. 
In the present case, the question would be the following: is the gravitational 
entropy computable from the Fourier coefficients of Maass-Hecke eigenformsl 
Another remark is the following. In the semi-classical limit, the expected 
statistics for the energy levels spacing distributions for X = PSL(2,Z)\M 
is the Poisson distribution. Indeed, from general considerations regarding 
chaotic systems, one would expect the system follow predictions from ran- 
dom matrix theory (GOE/GUE ensemble). These ensembles present level 
repulsion at short distance and rigidity at long distance. But, again, due 
to the arithmetic property, we have an anomalous statistics, the Poissonian 
one, which exhibits the clustering property. The knowledge of the expected 
statistics could, in principle, allow to compare with the observations. 
The inclusion of matter changes the shape of the billiard and also increases 
the number of dimensions of the billiard. The most interesting case is per- 
haps supergravity theory in 11 dimensions, which is believed to be a kind 
of ultimate theory unifying all fundamental interactions. The supergravity 
billiard is the fundamental Weyl chamber for E 10 , thus in this case the state- 
ment is that the wave function is a Maass waveform with respect to the 
(discrete) Weyl group of E 10 (which is not known). In this case, we can not 
speak safely of Maass cusp forms, and we must use only the general term 
Maass waveform, because we do not know if the residual spectrum is empty 
as in the cases of each T(N) (the congruence subgroups of PSL(2, Z)). Be- 
sides, the discrete spectrum could also be degenerate. The existence of Maass 
forms with respect to W(E 10 ) should be guaranteed by the fact W(E 10 ) is 
arithmetic. 

Usually in quantum mechanics, one has a discrete and a continuous spectrum 
for a self-adjoint Hamiltonian and these are separated. Bound states are 
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the proper eigenf unctions of the discrete spectrum, whereas the continuous 
part is interpreted as a free motion. Our model of quantum cosmology is 
different from this usual situation. In fact, from the Selberg theory for the 
automorphic Laplacian inside a fundamental domain of some T(N), we know 
that the discrete spectrum is not separated from the continuous one [7,00) 
(whose improper eigenfunctions are given by the Eisenstein series), but it is 
embedded in the continuous parfl We believe this says something about the 
nature of quantum gravity, i.e. this suggests that quantum gravity /cosmology 
is a non-trivial mixing of discrete and continuous concepts. 
It would be interesting to understand if other approaches to the problem of 
quantum gravity like the formalism of loop quantum gravity or string theory 
say something regarding this point. 

5 Scarring and Asymptotic Number of Quan- 
tum States 

In [TJ , Barrow and Levin have analyzed the case of a finite universe emerging 
from a compact octagon on the hyperbolic plane. The octagon is a com- 
pact domain (but the associated discrete group is not arithmetic), thus the 
classical motion inside is chaotic. They found a fractal structure, which, in 
the classical to quantum transition, can persist in forms of scars, ridges of 
enhanced amplitude in the semi-classical wave function. They conclude that 
if the universe is finite and negatively curved, the cobweb of luminous matter 
might be a residue of these primordial quantum scars. 

Of course we can never know if our universe is finite or not. Yet, in the case 
of a generic singularity, we have seen that the interesting domain is the one 
for PGL(2, Z), which is an arithmetic group. Thus, there is no scarring (see 
Appendix [B]). The conclusions of Barrow and Levin do not apply in this case 
and it remains to understand the physical (cosmological) meaning for the 
absence of these scarred states. Results in this direction by loop quantum 

11 There are cases in scattering problems where points of the discrete spectrum lie in 
the continuous spectrum, but our situation is distinguished, because the whole discrete 
spectrum is embedded in the continuous one. 
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gravity/cosmology or string theory would be interesting, in order to confirm 
or discard the roles of scarred stated in quantum cosmology 
Note one more thing. Our analysis is limited to the case of pure gravity 
in 4 dimensions, thus the conclusion about the absence of scarred states 
in quantum cosmology is valid only in this context. Thus one may think 
that the inclusion of matter or extra dimensions could change the situation. 
If we consider the case of supergravity theory in 11 dimensions (or a not 
well defined quantum version of it, like M-theory or whatever), then, as we 
mentioned, the cosmological billiard is the Weyl chamber of E 10 [6] (the role 
of the fermions is still matter of debate) whose Weyl group is still arithmetic. 
Thus for the hyperbolic manifold W / (Eio)\lH[ 9 the quantum unique ergodicity 
theorem should be true and, again, there is no scarring effect. 11-dimensional 
supergravity is a candidate theory to describe the universe and all of its 
interactions. The message is that, with the knowledge we have today, it does 
seem that in the early universe scarred states are absent. And, in our opinion, 
quantum gravity has to cope with quantum chaos. The absence of scarred 
states is deduced from the arithmetic quantum unique ergodicity theorem, 
which says that the only invariant measure at the quantum level is the one 
given by the Lebesgue measure on the hyperbolic manifold (which means 
uniqueness of the invariant measure for this model of quantum cosmology). 
In other words, all the features of classical chaos disappear at the quantum 
level, giving a system which is closer to a classically integrable system. 
Finally, regarding the asymptotic number of quantum states, we can count 
them using the Selberg result for PSL(2, Z) 



for R — > oo. This gives the asymptotic number of eigenvalues/eigenfunctions, 
thus it counts the asymptotic number of discrete states in quantum cosmology. 
Actually, we have to consider only the eigenvalues whose eigenf unctions are 
even, so only a part (on average, half) of this asymptotics is physical. Again, 
a confirmation of this from other approaches to the problem to quantum 
cosmology would be extremely interesting. For E 10 , there is no such a result. 
The first eigenvalues of A on X are E 1 = 91.12, E 2 = 148.43, E 3 = 190.13, E 4 = 




R 



(27) 



0<Ej<R 
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206.16. They are respectively even, odd, odd, even and odd, with respect to 
the symmetry about x — 0. See [12] and [20] for the description of the nu- 
merical methods used to determine eigenvalues and an approximate form for 
the eigenfunctions. It is interesting also to compare the energy levels com- 
puted in [18] using the approximate domain with the ones computed from 
the "exact" domain (but always via numerical methods). 

6 Conclusions, Some Speculations and Points 
to be developed in the future 

Assuming the billiard representation introduced by DHN in the study of 
the cosmological singularity in the BKL limit, we have commented on the 
question of integrability/chaoticity of the null billiard flow inside the Weyl 
chamber of HA^: one must be careful about statements which derive the 
chaoticity of cosmological billiard from the finite- volume argument (typical of 
the Euclidean case), as we are in a pseudo-Riemannian setting. In fact, from 
the Anosov property of the flow on the basis of the billiard (the fundamental 
domain for PGL(2, Z)), it does not follow the chaoticity for the dynamics in 
the full 3-dimensional Weyl chamber, since from the ergodicity of a dynamical 
system in a proper subset of the phase space one can not infer the ergodicity 
of the system in the full phase space. We do hope to study billiards and 
flows in pseudo-Riemannian manifolds in the near future, since these are the 
natural evolution of classical billiards towards considering the problem of 
ergodicity in special/general relativity. The next step, in fact, should be the 
study of relativistic billiards confined in some polyhedron in a Minkowski 
spacetime. 

We have also carried on the quantum analysis for the cosmological billiard 
using the arithmetic properties of the modular group. The result is that, for 
pure gravity in 4 dimensions, the wave function of the universe (or better, 
the angular part projected on the hyperbolic plane, i.e. the space part of the 
full wave function) is an automorphic form for PGL(2,Z), precisely an even 
Maass cusp form for PSL(2,Z). Indeed, since PSL(2,Z) is arithmetic, the 
wave function is a Maass-Hecke eigenform, being also eigenfunction of the 
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Hecke operators (which commute with the Hamiltonian) . In our opinion, the 
most important prediction/hint of the model is that the discrete spectrum 
is fully embedded in the continuous spectrum, which is a very distinguished 
situation. In fact, for every point in the discrete spectrum and every normal- 
izable even Maass cusp form, there exists a corresponding not-normalizable 
Eisenstein series with the same value of the generalized eigenvalue. This is 
an odd situation, thus we would like to understand if this phenomenon hides 
some deep truth which is an intrinsic property of quantum gravity/cosmology 
or if it is due to the billiard model which is too simplified. 
The arithmetic nature of PGL(2, Z) allows to state that in the early universe 
scarred states are absent, thanks to the quantum unique ergodicity theorem 
by Lindenstrauss. The model allows also to count the asymptotic number 
of quantum states, whose number is half of the full Selberg asymptotics 
for PSL(2, Z), and to identify a distribution for the level spacing statistics. 
This distribution is the Poissonian one and is characterized by the clustering 
property. One should understand the meaning of this as opposite to the level 
repulsion predicted from random matrix theory for non-arithmetic chaotic 
systems. 

These conclusions, true for pure gravity in 4 dimensions if we believe the bil- 
liard representation, are very likely valid also in the case of 11-dimensional 
supergravity close to the singularity, whose billiard is modelled on the E 10 
hyperbolic Kac-Moody algebra. Very briefly, one can say that the model 
predicts specific automorphic properties for the wave function, it gives the 
asymptotic number of quantum states and the statistics for the level spacing. 
Arithmetic systems at the quantum level seem to be closer to integrable sys- 
tems than to generic chaotic systems, thus, for this model, any manifestation 
of classical chaos disappears at the quantum level. 

In [7j, we have pointed out that the Selberg trace formula for PSL(2, Z) gives 
a semiclassical quantization rule for pure gravity in 4 dimensions, as typically 
occurs for the Gutzwiller trace formula in quantum chaos. But the difference 
is that the Gutzwiller trace formula is divergent, whereas the Selberg trave 
formula is convergent. Thus, in the BKL limit to the singularity and in the 
billiard representation, a semiclassical quantization of gravity is well defined. 
Indeed, the Selberg trace formula is a kind of path-integral. This point of 
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view is emphasized in the book by C. Grosche [Sj- But remember that indeed 
we have an infinite number of semi-classical quantization rules, since the class 
of test functions entering the trace formula is very large [20]. This is the best 
one can do in a rigorous/convergent way for the quantum systems whose 
semi-classical limit is a Hamiltonian flow for which the Selberg trace formula 
is valid. Nevertheless, the trace formula contains a lot of information, in fact 
Selberg derived his asymptotics using a particular Gaussian test function. 
It may be that this quantization helps in finding a correct theory of quan- 
tum gravity, as it happened in the early days of quantum mechanics after 
the introduction of Borh-Sommerfeld quantization rules. The Selberg trace 
formula, in a sense, is an identity between quantum mechanics (because it 
contains the spectrum of the quantum Hamiltonian) and classical mechanics 
(because it contains a sum over all classical periodic orbits). In [7], we have 
made an attempt to interpret the imaginary roots of the hyperbolic algebra 
in terms of the periodic orbits which enter the Selberg trace formula. These 
periodic orbits live inside PSL(2, Z), which is the projected billiard. The idea 
is to check if these periodic orbits on the hyperbolic plane can be lifted up 
to periodic orbits (which escape the singularity) in the 3-dimensional billiard 
such that they correspond to solutions (also approximate) of Einstein equa- 
tions in the asymptotic regime as it happens for the Kasner solutions. In 
fact, to the simple roots of the algebra and their Weyl-images (i.e. the real 
roots) one can associate Kasner solutions, because every reflection is a Weyl 
reflection from a Kasner epoch to another one. We would like to understand 
if something similar is possible for the imaginary roots, in particular relating 
the Selberg trace formula to the structure of the hyperbolic algebra. A future 
paper of more mathematical taste is in progress. 

The billiard representation is very interesting from many points of view. It is 
evident that in this formalism one has another example of a relation between 
general relativity and geometrical optics. One can wonder if this reformula- 
tion can be applied to other models too. We would like to understand if it 
is possible to reformulate problems of general relativity as billiard-like prob- 
lems (although billiard problems are generally very difficult to solve). This 
would help in studying the integrability or ergodicity of general relativity in 
particular situations. 
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For example, the variables used to describe Choptuik's solution to the col- 
lapse of a spherically symmetric scalar field show a profile with spikes which 
resemble (billiard) bounces, see for example figure 3 in [T0]H Moreover, this 
system exhibits discrete self-similarity, a property which is very likely shared 
by the BKL limit too (let us remember that the collapse studied numerically 
by Choptuik implies the formation of a black-hole, that is a singularity, thus 
possible analogies with the cosmological singularity are reasonable). How- 
ever, let us remember that Choptuik's solution is physical, while the BKL 
behavior, as far as we know, is theoretical. 

We have limited our analysis to the case of pure gravity in 3+1 dimensions, 
where, in the limit to the singularity, spatial points decouple and the dy- 
namics involves only a single (time) variable. The question now is: what 
could be the analog for a quantum theory when all 4 (time and space) vari- 
ables must be considered? Is it possible to describe gravity with the tools of 
higher-dimensional modular groups in particular regimes where one does not 
have spatial decoupling? Of course, this is a statement motivated more by 
beautiful mathematics than physics, thus it can be wrong. However, let us 
remember that higher- dimensional groups share with PSL(2, Z) a lot of nice 
properties, like existence of waveforms, trace formulae (the so-called Arthur 
trace formula) etc. In other words, the construction of a quantum theory 
would be viable under the the aegis of what we would call an "automorphic 
principle"!^. An interesting mathematical mechanism which could do the 
job is the lifting of automorphic forms to higher dimensions, but again the 

12 The author thanks prof. T Damour for explaining Choptuik's discover to him. 

13 Note that this proposal is different from the one contained for example in [6] (and 
references therein). DHN's approach is based on the assumption (for which there is indeed 
evidence) that supergravity theory in 11 dimensions (or a quantum extension of it) exhibits 
a hidden Eio symmetry at the level of the Lagrangian, i.e. there should exist a formulation 
of this theory invariant under the infinite-dimensional Lie group Eio\K(Eio), where K(Eio) 
is the formal maximal compact subgroup of the Eio Lie group. In this setting, the discrete 
group one looks for to build the quantum theory is a discretized version of Eio, Eio(Z). Our 
proposal is different: we suggest to look at discrete groups like GL(n, Z) and subgroups, 
not to some HA^(Z). Building automorphic forms for GL(n,Z) groups (and subgroups) 
is a well-established technique, while automorphic forms for Eio(Z) or HA^(Z), at the 
moment, are out of reach. Anyway, all these considerations can have nothing to do with 
true physics. 
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connection with the physics is all to prove. Everything is mathematically 
well defined and beautiful, and leads, unavoidably, to the classical Langlands 
program. Very recently, E. Witten et al have given a physical interpreta- 
tion of the geometric Langlands program in terms of gauge theory. This 
construction uses a lot Hecke eigensheaves, whose classical counterpart is 
Hecke operators and Hecke eigenfunctions. It is remarkable, in my opinion, 
that these objects appear in general relativity in a different context and with 
a different language. It could be another subtle indication of the so-called 
gauge-gravity correspondence, realized via primitive (arithmetic) objects. 

7 Acknowledgements 

The author wishes to thank profs. T. Damour, F. Englert, L. Houart, M. 
Henneaux for fruitful discussions at early stages of this work, prof. G. Espos- 
ito for reading a copy of the manuscript and prof. A. Sciarrino for suggestions 
and collaboration in the years. Finally, the author wishes to thank the ref- 
erees for their useful comments and for pointing out imprecisions and/or 
improvements. 



A The Hyperbolic Kac- Moody Algebra HA 



In this section we do not review the theory of Kac-Moody algebras (see [H] 
for details), but just fix the notation. 

The algebra HA^ is an infinite-dimensional Lie algebra, being a hyperbolic 
Kac-Moody algebra. This means that it is identified by the following (gen- 
eralized) Cartan matrix 




(28) 



or equivalently by its Dynkin diagram 
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The precise definition of the algebra is standard and goes as follows. Let f) 
be a complex vector space whose dimension is 3, and f)* its dual. Then there 
exist linearly independent indexed sets 

n := { ai } C f)* , and n v := {^} C f) , (29) 

such that aj (hi) = (i,j = 1,2,3). The Oj (hi) are called simple roots 
(dual simple roots). The sets II and II V are uniquely determined by A up to 
isomorphism. 

Then the Kac-Moody algebra HA^ is the complex Lie algebra generated by 
f) U {ej, fi} with the following defining relations: 

ft] = 5ij hi , [h, h'\ = for h,h'ef) (30) 

\hii Cj] Oiij 6j , \hi) fj\ Q>ij fi 

(ad ei y- a - ej = , (ad /O 1 "^ = for ^ j . 

[) is the maximal abelian subalgebra of q(A) (Cartan subalgebra). The order 
of the Cartan matrix is also called the rank of the algebra, which in this case 
coincides with the rank of the matrix A. Finally, the algebra is simple. 
The term hyperbolic refers to the fact that its Cartan matrix has determinant 
— 1 and Lorentzian signature ( — h+). In a sense, it is the simplest hyperbolic 
Kac-Moody algebra of rank 3, being the canonical hyperbolic extension of 
Ai (su(2)) trough the affme node a 2 and the hyperbolic node a 3 . The scalar 
products between the simple roots ( 

(ai,a 2 ) = -2, (0:1,0:3) =0, (o 2 ,o 3 ) = -l (31) 

A is a flat metric on a certain pseudo-Riemannian vector space fjj^ (see be- 
low). Given the simple roots, one has the notion of real roots, which are 
W— equivalent to the simple roots, i.e. they are the image through the Weyl 
group of the simple roots, and imaginary roots. It turns out that the real 
roots a have positive squared norm, while for imaginary roots a 2 < 0. In- 
deed, for hyperbolic (symmetrizable) Kac-Moody algebras, the imaginary 
roots are precisely all the vectors in the root lattice ® i Z«j which have zero 
or negative squared norm (Moody theorem). 
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The Weyl group W is defined by the fundamental reflections in the simple 
roots 

R i (P)=P-2p^a t (32) 

It is a Coxeter group, i.e. it is a discrete group generated by the fundamental 
reflections with the following relations 

R\ = 1 , (RiRjT 13 = 1 (33) 

where the Coxeter exponents are positive integers or oo (in this case 
one puts x°° = 1). For Kac-Moody algebras, one has = 2,3,4,6 or oo 
according to a^- = 0, 1, 2, 3 or > 4. Thus, different Kac-Moody algebras 
with the same Coxeter exponents has isomorphic Weyl group M For HAf 

R 2 t = (RM 2 = (R 2 R 3 ) 3 = 1 (35) 

and the Weyl group can be easily identified with PGL(2, Z) putting T = 
R2R1 and S = R1R3, where T, S are the standard generators for PSL(2,Z). 
The even (positive) Weyl group W + is thus isomorphic to PSL(2, Z). The 
properties of PSL(2, Z), the modular group, are recalled in the next section. 
Every Coxeter group defines a Coxeter polytope, which is the polyhedron 
whose faces are point-wise fixed by the fundamental reflections. When the 
Coxeter group is the Weyl group of a Kac-Moody algebra, this polyhedron 
is called the Weyl chambeil 15 ! of the algebra and is contained in fjjjj . Thus 
the faces of the Weyl chamber are orthogonal (with respect to the metric 
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For example, the two rank-3 algebras 





(34) 



have the same Weyl group as HA^. 

15 To be honest, the Weyl chamber is the domain in f) defined by 

{he t) R \ai{h) >0} (36) 

where fjjg. = {h € f)| a, (ft) G 8} is a W— stable real subspace of rj. In this notation, what 
we call Weyl chamber should more correctly called the dual Weyl chamber, since we are 
interested in the polyhedron in f)Jj. 
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given by the Cartan matrix) to the simple roots and are hyperplanes in the 
pseudo-Riemannian space f)Jj. 

B Maass Automorphic Forms for PSL(2, Z) 

Let us remind the main properties of the modular group PSL(2, Z). First we 
consider PGL(2, Z), sometimes called the extended modular group, whose 
fundamental domain we call JF 3 (remember that one of the angles is 7r/3). It 
is a discrete group acting on the hyperbolic plane and is generated by the 
hyperbolic reflections in the three sides 

R 1 (z) = -z, R 2 (z) = -z+l, i? 3 (z) = = (37) 

PSL(2,Z) is its most important subgroup, its domain T is twice the funda- 
mental domain of PGL(2, Z) 

T= jzeH: |*| > l,|x| < ij (38) 
and its hyperbolic area is 

ti?) = f ( 39 ) 

where \x is the usual hyperbolic measure on the hyperbolic plane, d\i = 
dxdy/y 2 . The standard generators for PSL(2,Z) are 

considered as fractional linear transformations, i.e. 

T{z) =z + l, S(z) = -- (41) 

z 

S is such that S 2 = 1 (S = S' 1 ); moreover Tr S = < 2, Tr T = 2 thus 
S is an elliptic element, T is a parabolic element (for details on hyperbolic 
geometry and Fuchsian groups see [15]). 

T and S can be expressed of course in terms of the reflections Ri belonging 
to PGL(2,Z) 

T = R 2 R 1 ^ RiR 2 = T~\ S = R\R 3 = R3R1 (42) 
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Figure 4: The modular domain (the domain for PGL(2, Z) is the desym- 
metrized one) and the modular surface (from [20J). 

Note that PGL(2,Z) contains hyperbolic reflections, that is negative isome- 
tries with determinant —1, while PSL(2, Z) contains only positive isometries. 
Thus we can consider the hyperbolic surface X := PSL(2, Z)\H, called the 
modular surface, which is an oriented, finite area, non-compact 2-dimensional 
Riemannian manifold. It is also a Riemann surface with the complex struc- 
ture inherited by H. It has a cusp at infinity, which corresponds to the fixed 
point ioo of the parabolic transformation T. We are interested in the spectral 
problem for X. See the excellent survey (20] for more details and references. 
The so-called fundamental spectral problem of quantum chaos is the following 



Here A = y 2 [d 2 x + d 2 ) is the hyperbolic Laplacian. The numbers = Eq < 
Ei < E 2 < • • • for the which the spectral problem has solutions form the 
discrete spectrum (eigenvalues) for X. The only known eigenvalue is E = 0, 
to which it corresponds the constant eigenfunctions ip { z ) (remember that 
T has finite area). Solutions to the previous spectral problem are called 
Maass automorphic forms (or Maass waveforms), after the mathematician 
H. Maass, who first studied them. The term automorphic refers to the peri- 
odicity under PSL(2,Z), i.e. tpijz) = ip(z). It is not evident that solutions 
to this problem exist, since the modular surface is not compact. Moreover, 
no explicit eigenvalues or eigenfunctions are known or expected for X, al- 




(43) 
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though Maass himself produced an explicit subsequence of eigenvalues for 
other discrete groups (e.g. T(4)). 

One would like to formulate for this case the analog of the Weyl law for the 
Dirichlet problem for the Euclidean Laplacian on a compact domain in M 2 
on the asymptotic number of eigenvalues. This is possible thanks to Selberg 
theory, in particular his famous trace formula. Selberg found that there 
exists a continuous spectrum too, and the latter is the whole interval [| , oo) 
with multiplicity one. The corresponding not-normalizable eigenfunctions 
are given by the Eisenstein series, which, for X, read as follows 



E(z,s)= forRe S >l (44) 



where is the stabilizer of oo, which is the only cusp, z = x + iy and the 
transformation g is represented as usual by the matrix ( \ J . The Eisen- 



c d 

stein series extend meromorphically to the whole C— plane and are analytic 
on Re s = \. The continuous spectrum is thus furnished by these generalized 
eigenfunctions E (z, \ + it) , t > 0, 

A£^,i + ^ + Q + /) E^ + itj =0 (45) 

and these are of course PSL(2, Z)— periodic 

E(gz, s) = E(z, s) for any g £ PSL(2, Z) (46) 

that is they are generalized Maass automorphic forms. Let us call ip(s) 
the constant term in the Fourier expansion of the Eisenstein series: this is 
meromorphic in C and for Re s > 1/2 its poles are in (|, l] . The residues 
at these poles are solutions to the spectral problem and form the residual 
spectrum of X. If we now take the orthogonal complement in L 2 (X,fi) of 
the continuous and residual spectrum, we obtain the cuspidal space Ll ups (X). 
It is invariant under the hyperbolic Laplacian and the resolvent (A — A) -1 
is compact when restricted Ll usp . A Maass form which also lies in L 2 cusp is 
called a Maass cusp form. These cusp forms are a kind of building blocks 
for the theory of automorphic forms. Their existence is of course tied to the 
size of L 2 cusp , since L 2 cusp ^ {0} is not obvious. 
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However, it turns out that for PSL(2, Z), ip(s) has no poles in f~, l), thus the 
residual spectrum is empty and any Maass form is automatically a cusp form. 
Using his trace formula and the fact that for PSL(2, Z) ip(s) can be expressed 
in terms of a ratio involving the (completed) Riemann zeta-function, Selberg 
was able to show that the contribution of the continuous spectrum to the 
Weyl law is negligible, and that 

0<Ej<R 

that is solutions to the spectral problem exist and in abundance (each eigen- 
value is counted with its multiplicity). Let us stress that no eigenvalues/eigenfunctions 
are explicitly known. 

The existence of Maass cusp forms for PSL(2,Z) is deeply related to the 
its arithmetic nature^. One can also consider the spectral problem for 
PGL(2,Z); this latter contains negative isometries, thus one can not form 
a hyperbolic manifold by taking the quotient, but one can consider the spec- 
tral problem with Neumann boundary conditions 

- Aip = Eip 
ip e L 2 {F 3 ,fi) 

d n ipW s = (48) 

Solutions to this problem are given by the even Maass cusp forms for PSL(2, Z). 
In full generality, let X(N) be the surfaces built out of the congruence sub- 
groups r(iV). The spectral problem applies exactly as above. Regarding the 
low-energy spectrum of X(N), let us call Ei(N) the closest eigenvalue to 
E = 0. A deep conjecture (still open) of Selberg states that 

Ei(N) > \ (49) 

For N = 1 (X = X(l)), the case of the modular group, the numerical 
calculations show that the first eigenvalue is E\ = 91.12 . . .. Since there is no 

16 Phillips and Sarnak have more generally shown that the arithmetic groups are the 
only ones which allow for solutions. Many numerical experiments support the absence 
of eigenvalues and Maass forms for non- arithmetic cases. For the precise definition of an 
arithmetic Fuchsian group see [15] . 
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residual spectrum! 17 !, we can take for Ei the smallest eigenvalue of a Maass 
cusp form on X. One can understand the number 1/4 by recalling that a 
result due to McKean shows that the spectrum of A on the universal covering 



of A on X is simple (the situation may indeed be very different for the other 
surfaces X(N)), a result confirmed by many numerical experiments. For 
the high-energy behavior of the spectrum, one would expect, using standard 
arguments from random matrix theory, that it would fit the GOE-ensemble 
predictions (the system is time- invariant), as the geodesic flow on X is an 
Anosov flow (thus chaotic). Yet, all the numerical experiments show that 
the high-energy spectrum follows a Poissonian distribution (which is instead 
typical of systems whose semi-classical limit is integrable). The reason for 
that is, again, the arithmetic nature of PSL(2, Z), in particular the existence 
of additional symmetries, the so-called Hecke operators, which for n > read 
as follows 



the sum going over all positive integers a, d with ad = n and < b < d. These 
operators are a commutative familyf^l of self-adjoint operators on L 2 (X) and 
commute with A and the reflection R\. Thus they preserve the even/odd 
eigenspaces of A and each eigenspace has a basis consisting of simultaneous 
eigenfunctions of all Hecke operators. Such eigenfunctions are called Maass- 
Hecke eigenforras. Given such an eigenfunction if>, T n ip = \ n ip, its Fourier 
coefficients are given by 



and we can normalize the first Fourier coefficient a^(l) = 1; in this way the 
n— th Fourier coefficient is the Hecke eigenvalue A n . Hecke eigenvalues enjoy 

17 besides E = 0, which thus should not be considered in the discrete spectrum. 

18 Mathematicians have no doubts that the anomalous statistics for arithmetical dynam- 
ical systems is due to these Hecke operators, which at the quantum level commute with 
the Laplacian A and somehow mimic an integrable system. A rigorous proof of that is 
lacking. 



L 2 (H) is [J, oo). Moreover, P. 



Cartier conjectured that the cuspidal spectrum 




(50) 




(51) 
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a lot of nice properties, the most important is that they are multiplicative 




A 



•rnn 



X m \ n m, n co-prime 
Apfc+i + X p k-i p prime 



(52) 
(53) 



Let us conclude this section by recalling a deep and recent mathematical 
result about the quantum unique ergodicity for arithmetic hyperbolic mani- 
folds. Let us first observe that, from the classical point of view, arithmetical 
dynamical systems are chaotic as any other model on compact negatively 
curved manifolds. But at the quantum level the arithmetic property is pe- 
culiar. If one defines the probability measures /ij 



where ipj is an eigenfunction of the Laplacian and dvol is the Riemannian vol- 
ume element, then the quantum unique ergodicity conjecture by Rudnick and 
Sarnak states that for compact manifolds of negative curvature, the measures 
fij converge to dvol (in the weak* topology). The conjecture by Rudnick and 
Sarnak, if true, is remarkable, because it asserts that at quantum level and in 
the semi-classical limit, there is no manifestation of chaos. In particular, one 
would have quantum unique ergodicity, that is only one possible quantum 
limit, whereas classical unique ergodicity, i.e. uniqueness of the invariant 
measure for the Hamiltonian flow, is never satisfied for chaotic systems. 
This conjecture is motivated by a theorem due to Shnirelman, which says that 
for an ergodic system there exists a subsequence jk for which the measures 
Hj k converge to the standard normalized Lebesgue measure. Shnirelman's 
theorem is an statement of equidistribution, that is one often says that the 
eigenfunctions equidistribute because the probability densities \ipj k \ tend to 
a constant independent on any point on the manifold. Thus, for non excep- 
tional E n , \ip n \ 2 can never localize to just a finite number of closed geodesies. 
This theorem has been improved by Colin de Verdiere and Zelditch, who 
showed that eigenfunctions still equidistribute on non-compact manifolds like 
PSL(2,Z)\H. The presence of an exceptional set is of course troubling. In 
stadium-like domains, the so-called scarring effect has been observed. For 
numerous n, the topography of ip n is found to contain clear ridges of mass 
of scars, situated roughly along what would appear of closed geodesies. The 



d[ij = \ipj\ 2 dvo\ 



(54) 
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location of these scats changes with n. In practice, scars is what is left of 
periodic orbits [12] . 

But for negatively curved manifolds, the QUE conjecture by Rudnick and 
Sarnak denies the existence of these scars. The proof of this conjecture is still 
lacking, yet progress has been made for arithmetic hyperbolic surfaces thanks 
to E. Lindenstrauss. He has shown that for a compact arithmetic quotient the 
quantum unique ergodicity conjecture is true and that, moreover, the same 
statement holds for the modular surface (which is non-compact). Thus, in 
the latter case, there is no scarring, i.e. these is no manifestation of classical 
chaos at the quantum level. 
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